The QGP dynamics in relativistic heavy-ion coUisions 

. E. L. Bratkovskaya^'^, W. Gassing'^, V. P. Konchakovski^'^ , 

O. Linnyk^, V. Ozvenchuk^, V.D. Toneev^'^ and V. Voronyuk^'^'^ 

c 2 ^ ' ^ Institute for Theoretical Physics, University of Frankfurt, Frankfurt, Germany 

Cn , ■^ Frankfurt Institute for Advanced Study, Frankfurt, Germany 

$H ' ^ Institute for Theoretical Physics, University of Giessen, Giessen, Germany 

C*h' * Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 

■^r , ^ Joint Institute for Nuclear Research, Dubna, Russia 

m 

Abstract. The dynamics of partons and hadrons in relativistic nucleus- nucleus collisions is 

analyzed within the novel Parton-Hadron-String Dynamics (PHSD) transport approach, which is 

based on a dynamical quasiparticle model for the partonic phase (DQPM) including a dynamical 

hadronization scheme. The PHSD model reproduces a large variety of observables from SPS 

to LHC energies, e.g. the quark-number scaling of elliptic flow, transverse mass and rapidity 

O ■ spectra of charged hadrons, dilepton spectra, open and hidden charm production, collective flow 

2 ' coefficients etc., which are associated with the observation of a sQGP. The 'highlights' of the 

, ^ , , latest results on collective flow are presented and open questions/perspectives are discussed. 

> 

1. Introduction 

The dynamics of the early universe in terms of the 'Big Bang' may be studied experimentally by 
ultrarelativistic nucleus-nucleus collisions at Relativistic-Heavy-Ion-CoUider (RHIC) or Large- 
f~^ I Hadron-Collider (LHC) energies in terms of 'tiny bangs' in the laboratory. The Power Spectrum 

C^ ' extracted from the Cosmic Microwave Background Radiation has some analogy to the Fourier 

components of particles in the azimuthal angular distribution \T\. The discovery of large 
azimuthal anisotropic flow at RHIC has provided conclusive evidence for the creation of dense 
jv>( ■ partonic matter in ultra-relativistic nucleus-nucleus collisions. With sufficiently strong parton 

^ ', interactions, the medium in the collision zone can be expected to achieve local equilibrium and 
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exhibit approximately hydrodynamic flow O El H]. The momentum anisotropy is generated 
due to pressure gradients of the initial "almond-shaped" collision zone produced in noncentral 
collisions [21 [5] . The azimuthal pressure gradient extinguishes itself soon after the start of the 
hydrodynamic evolution, so the final flow is only weakly sensitive to later stages of the fireball 
evolution. The pressure gradients have to be large enough to translate an early asymmetry 
in density of the initial state to a final-state momentum-space anisotropy. In these collisions 
a new state of strongly interacting matter is created, being characterized by a very low shear 
viscosity r] to entropy density s ratio, r]/s, close to a nearly perfect fluid [5l[6l[7j. Lattice QCD 
(IQCD) calculations [8, 9j indicate that a crossover region between hadron and quark-gluon 
matter should have been reached in these experiments. 

An experimental manifestation of this collective flow is the anisotropic emission of charged 
particles in the plane transverse to the beam direction. This anisotropy is described by the 
different flow parameters defined as the proper Fourier coefficients v^ of the particle distributions 
in azimuthal angle ip with respect to the reaction plane angle ^rp- At the highest RHIC collision 



energy of y/s^N = 200 GeV, differential elliptic flow measurements V2{pt) have been reported 
for a broad range of centralities or number of participants Npart- For Npart estimates, the 
geometric fluctuations associated with the positions of the nucleons in the collision zone serve as 
the underlying origin of the initial eccentricity fluctuations. These data are found to be in accord 
with model calculations that an essentially locally equilibrated quark gluon plasma (QGP) has 
little or no viscosity |10^ [TTl [T2] . Collective flow continues to play a central role in characterizing 
the transport properties of the strongly interacting matter produced in heavy-ion collisions at 
RHIC and even LHC and shed some light on the scale of initial state fluctuations. 

The Beam-Energy-Scan (BES) program proposed at RHIC |13j covers the energy interval 
from y/SNN = 200 GeV, where partonic degrees of freedom play a decisive role, down to the 
AGS energy of s/snn ~ 5 GeV, where most experimental data may be described successfully in 
terms of hadronic degrees-of- freedom, only. Lowering the RHIC collision energy and studying 
the energy dependence of anisotropic flow allows to search for the possible onset of the transition 
to a phase with partonic degrees-of-freedom at an early stage of the collision as well as possibly 
to identify the location of the critical end-point that terminates the cross-over transition 
at small quark-chemical potential to a first order phase transition at higher quark-chemical 
potential [HlfTS] , 

This contribution aims to summarize excitation functions for different harmonics of the 
charged particle anisotropy in the azimuthal angle at midrapidity in a wide transient energy 
range, i.e. from the AGS to the top RHIC energy. The first attempts to explain the preliminary 
STAR data with respect to the observed increase of the elliptic flow V2 with the collision energy 
have failed since the traditional available models did not allow to clarify the role of the partonic 
phase |16j . In this contribution we investigate the energy behavior of different flow coefficients, 
their scaling properties and differential distributions (cf. Ref. [171 E]). Our analysis of the 
STAR/PHENIX RHIC data - based on recent results of the BES program - will be performed 
within the Parton-Hadron-String Dynamics (PHSD) transport model [19] that includes explicit 
partonic degrees-of-freedom as well as a dynamical hadronization scheme for the transition from 
partonic to hadronic degrees-of-freedom and vice versa. For more detailed descriptions of PHSD 
and its ingredients we refer the reader to Refs. ^0\ EU [221 123]. 

2. Results for collective flows 

We directly continue with the results from PHSD in comparison with other approaches and the 
available experimental data. 

2.1. Elliptic flow 

The largest component, known as elliptic flow V2, is one of the early observations at RHIC [24J. 
The elliptic flow coefficient is a widely used quantity characterizing the azimuthal anisotropy of 
emitted particles, 

2 2 

V2 =< cos{2i; - 2*) >=< ^;~^^ > , (1) 

Px + Py 

where ^ rp is the azimuth of the reaction plane, px and py are the x and y component of the 
particle momenta and the brackets denote averaging over particles and events. This coefficient 
can be considered as a function of centrality, pseudorapidity rj and/or transverse momentum px. 
We note that the reaction plane in PHSD is given by the (x — z) plane with the z-axis in beam 
direction. The reaction plane is defined as a plane containing the beam axes and the impact 
parameter vector. 

We recall that at high bombarding energies the longitudinal size of the Lorentz contracted 
nuclei becomes negligible compared to its transverse size. The forward shadowing effect then 
becomes negligible and the elliptic flow fully develops in-plane, leading to a positive value of 
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Figure 1. (l.h.s.) The average elliptic flow V2 of charged particles at midrapidity for minimum 
bias collisions at ^/SNN = 9.2, 19.6, 62.4 and 200 GeV (stars) is taken from the data compilation 
of Ref. |16]). The corresponding results from different models are compared to the data and 
explained in more detail in the text, (r.h.s.) Evolution of the parton fraction of the total energy 
density at midrapidity (from PHSD) for different collision energies at impact parameters b =1 
fm and 10 fm. 



the average flow V2 since no shadowing from spectators takes place. In Fig. [T] (l.h.s.) the 
experimental V2 data compilation for the transient energy range is compared to the results from 
various models: PHSD, HSD as well as from UrQMD amd AMPT as included in Ref. [16j . The 
centrality selection is the same for the data and the various models. 

The HSD [271 [28l [29] and UrQMD (Ultra relativistic Quantum Molecular Dynamics) [251 126] 
are the hadron-string models and, thus, essentially provide information on the contribution 
from the hadronic phase |30j . As seen in Fig. [H being in agreement with data at the lowest 
energy y^J/vJv = 9.2 GeV, the HSD and UrQMD model results then either remain approximately 
constant or decrease slightly with increasing ^/sNN and do not reproduce the rise of V2 with the 
collision energy as seen experimentally. 

The AMPT (A Multi Phase Transport model) [31ll32] uses initial conditions of a perturbative 
QCD (pQCD) inspired model which produces multiple minijet partons according to the number 
of binary initial nucleon-nucleon collisions. The string melting (SM) version of the AMPT model 
(labeled in Fig. [1] as AMPT-SM) is based on the idea of 'melting' of hadrons or strings above 
the critical energy density of e ~ 1 GeV/fm^ to massless partons. The subsequent scattering 
of the quarks are based on a parton cascade with (adjustable) effective cross sections which are 
significantly larger than those from pQCD [31[ [32 j . Once the partonic interactions terminate, 
the partons hadronize through the mechanism of parton coalescence. 

We find from Fig. [1] that the interactions between the minijet partons in the AMPT model 
indeed increase the elliptic fiow significantly as compared to the hadronic models UrQMD and 
HSD. An additional inclusion of interactions between partons in the AMPT-SM model gives rise 
to another 20% of V2 bringing it into agreement (for AMPT-SM) with the data at the maximal 
collision energy. So, both versions of the AMPT model indicate the importance of partonic 
contributions to the observed elliptic flow V2 but do not reproduce its growth with y^J/vTv- The 
authors address this result to the partonic-equation-of state (EoS) employed which corresponds 
to a massless and noninteracting relativistic gas of particles. This EoS deviates severely from 
the results of lattice QCD calculations for temperatures below 2-3 Tc- Accordingly, the degrees- 
of-freedom are propagated without self-energies and a parton spectral function. 
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Figure 2. Average anisotropic flows v^ (l.h.s.) and V4 (r.h.s.) of charged particles at mid- 
pseudorapidity for minimum bias Au + Au collisions calculated within the PHSD (solid lines, 
red) and HSD (dashed lines, blue) models. 



The PHSD approach incorporates the latter medium effects in line with a IQCD equation-of- 
state and also includes a dynamical hadronization scheme based on covariant transition rates. 
As has been demonstrated in Refs. [171 [T8] and explicitly shown in Fig. [1] (l.h.s.), the elliptic 
flow V2 from PHSD (red line) agrees with the data from the STAR collaboration and clearly 
shows an increase with bombarding energy. 

An explanation for the increase in V2 with collision energy is provided in Fig. [1] (r.h.s.) where 
the partonic fraction of the energy density is shown with respect to the total energy where the 
energy densities are calculated at mid-rapidity. As discussed above the main contribution to 
the elliptic flow is coming from an initial partonic stage at high y/s. The fusion of partons to 
hadrons or, inversely, the melting of hadrons to partonic quasiparticles occurs when the local 
energy density is about e ~ 0.5 GeV/fm'^. As follows from Fig. [U the parton fraction of the 
total energy goes down substantially with decreasing bombarding energy while the duration of 
the partonic phase is roughly the same. The maximal fraction reached is the same in central and 
peripheral collisions but the parton evolution time is shorter in peripheral collisions. One should 
recall again the important role of the repulsive mean-field for partons in the PHSD model that 
leads to an increase of the flow V2 with respect to HSD predictions (cf. also Ref. [33j)- We point 
out in addition that the increase of V2 in PHSD relative to HSD is also partly due to the higher 
interaction rates in the partonic medium because of a lower ratio of rj/s for partonic degrees- 
of-freedom at energy densities above the critical energy density than for hadronic media below 
the critical energy density 134 ^ 135 ^ 136]. The relative increase in ^3 and V4 in PHSD essentially is 
due to the higher partonic interaction rate and thus to a lower ratio r]/s in the partonic medium 
which is mandatory to convert initial spacial anisotropies to final anisotropics in momentum 
space |3j 



2.2. Higher- order flow harmonics 

Depending on the location of the participant nucleons in the nucleus at the time of the collision, 
the actual shape of the overlap area may vary: the orientation and eccentricity of the ellipse 
defined by the participants fluctuates from event to event. Note, however, that by averaging 
over many events an almond shape is regained for the same impact parameter. 

Recent studies suggest that fluctuations in the initial state geometry can generate higher- 
order flow components [H \IU[ [55} I39j . The azimuthal momentum distribution of the emitted 
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Figure 3. (l.h.s.) Impact parameter dependence of anisotropic flows of charged particles at 
mid-pseudorapidity for minimum bias collisions of Au+Au at y/sNN = 200 GeV. Experimental 
points are from Ref. [33]. (r.h.s.) Beam energy dependence of the ratio Vi/{v2)^ for Au+Au 
collisions. The solid and dashed curves are calculated within the PHSD and HSD models, 
respectively. 



particles is commonly expressed in the form of a Fourier series as 
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where Vn is the magnitude of the n-th order harmonic term relative to the angle of the initial- 
state spatial plane of symmetry ^„. The anisotropy in the azimuthal angle tp is usually 
characterized by the even order Fourier coefficients with the reaction plane ^n = ^rp'- Vn = 
{exp{in{'tp — ^Rp))) [n = 2, 4, ...), since for a smooth angular profile the odd harmonics vanish. 
For the odd components, e.g. ^3, one should take into account event-by-event fluctuations with 
respect to the participant plane ^„ = ^ pp. We calculate the v^, coefficients with respect to ^3 
as: ?;3{^3} = (cos(3[^ — ^3]))/i?es(^3). The event plane angle ^3 and its resolution Res{^2,) 
are calculated as described in Ref. ^^ via the two-sub-events method f41'. 142]. 

In Fig. [2] we display the PHSD and HSD results for the anisotropic fiows U3 and ^4 of charged 
particles at mid-pseudorapidity for Au+Au collisions as a function of y/sJm. The pure hadronic 
model HSD gives U3 ~ for all energies. Accordingly, the results from PHSD (dashed red line) 
are systematically larger than from HSD (dashed blue line). Unfortunately, our statistics are 
not good enough to allow for more precise conclusions. The hexadecupole flow V4, stays almost 
constant in the energy range ^JsMN ^ 10 GeV; at the same time the PHSD gives noticeably 
higher values than HSD which we attribute to the higher interaction rate in the partonic phase. 

Alongside with the integrated flow coefficients Vn the PHSD model reasonably describes their 
distribution over centrality or impact parameter h. A specific comparison at s/snn = 200 GeV is 
shown in Fig. [3] for V2,V3 and V4. While V2 increases strongly with b up to peripheral collisions, 
V3 and f4 are only weakly sensitive to the impact parameter. The triangular fiow is always 
somewhat higher than the hexadecupole flow in the whole range of impact parameters b. 



2. 3. Ratios of different harmonics 

Different harmonics can be related to each other. In particular, hydrodynamics predicts that 
V4 oc (^2)^ [II]- The simplest prediction that V4 = 0.5(f2)^ is given for a boosted thermal freeze- 
out distribution of an ideal fluid, Ref. |45j . In this work it was noted also that ^4 is largely 
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Figure 4. (l.h.s.) Participant number dependence of the ^4/(^2)^ ratio of charged particles for 
Au+Au {s/snn = 200 GeV) colhsions. The experimental data points for 0.5< pT <0.8 GeV/c 
are from Ref. |4_6 j . (r.h.s.) Transverse momentum dependence of the ratio ^4/(1^2)^ of charged 
particles for Au+Au (at ^/snn = 200 GeV) collisions. The dashed and dot-dashed lines are 
calculated within the hydrodynamic approaches from Refs. [H] and [IS], respectively. The 
shaded region corresponds to the results from the AMPT model j48j . The experimental data 
points are from the STAR Collaboration |47j. 



generated by an intrinsic elliptic flow (at least at high pt) rather than the fourth order moment 
of the fluid flow. This is a motivation for studying the ratio ^4/(^2)^ rather than ^4 alone. As is 
seen in Fig. 4 (r.h.s.), indeed the ratio calculated within the PHSD model is practically constant 
in the whole range of ^/s^n considered but significantly deviates from the ideal fluid estimate 
of 0.5. In contrast, neglecting dynamical quark-gluon degrees-of-freedom in the HSD model, we 
obtain a monotonous growth of this ratio. 

The dependence of the ^14/(^2)^ ratio versus the number of participants Npart is shown in 
Fig. m for charged particles produced in Au + Au collisions at ^/sNN = 200 GeV. The PHSD 
results are roughly in agreement with the experimental data points from Ref. [47J but overshoot 
them for Npart ~ 250. 

As pointed out before, the ratio ^4/(^2)^ is sensitive to the microscopic dynamics. In this 
respect we show the transverse momentum dependence of the ratio ^4/(^2)^ in Fig. S] for 
charged particles produced in Au+Au collisions at y^s/wv = 200 GeV (20-30% centrality). 
The PHSD results are quite close to the experimental data points from Ref. |T7], however, 
overestimate the measurements by up to 20%. The hydrodynamic results - plotted in the same 
figure - significantly underestimate the experimental data and noticeably depend on viscosity. 
The partonic AMPT model j48] discussed above also predicts a slightly lower ratio than the 
measured one, however, being in agreement with both hydrodynamic models for p^ 0.8 GeV/c. 
Our interpretation of Fig. [J] (r.h.s.) is as follows: the data are not compatible with ideal 
hydrodynamics and a finite shear viscosity is mandatory (in viscous hydrodynamics) to come 
closer the experimental observations. The kinetic approaches AMPT and PHSD perform better 
but either overestimate (in AMPT) or slightly underestimate the scattering rate of soft particles 
(in PHSD). An explicit study of the centrality dependence of these ratios should provide further 
valuable information. 



3. Conclusions 

In summary, relativistic collisions of Au+Au from y^iivTV = 5 to 200 GeV have been studied 
within the PHSD approach which includes the dynamics of explicit partonic degrees-of-freedom 



as well as dynamical local transition rates from partons to hadrons and also the final hadronic 
scatterings. Whereas earlier studies have been carried out for longitudinal rapidity distributions 
of various hadrons, their transverse mass spectra and the elliptic flow V2 as compared to available 
data at SPS and RHIC energies \\.9\ 120) . here we have focussed on the PHSD results for 
the collective flow coefficients V2 , v^ and U4 in comparison to experimental data in the large 
energy range from the RHIC Beam-Energy-Scan (BES) program as well as different theoretical 
approaches ranging from hadronic transport models to ideal and viscous hydrodynamics. We 
mention explicitly that the PHSD model from Ref. |2U] has been used for all calculations 
performed in this study and no tuning (or change) of model parameters has been performed. 

We have found that the anisotropic flows - elliptic V2, triangular v^, hexadecapole V4 ~ are 
reasonably described within the PHSD model in the whole transient energy range naturally 
connecting the hadronic processes at lower energies with ultrarelativistic collisions where the 
quark-gluon degrees of freedom become dominant. The smooth growth of the elliptic flow V2 
with the collision energy demonstrates the increasing importance of partonic degrees of freedom. 
Other signatures of the transverse collective flow, the higher-order harmonics of the transverse 
anisotropy ^3 and ^4 change only weakly from ^/snn ~ 7 GeV to the top RHIC energy of 
^/snn = 200 GeV, roughly in agreement with experiment. As shown in this study, this success 
is related to a consistent treatment of the interacting partonic phase in PHSD whose fraction 
increases with the collision energy. 

The analysis of correlations between particles emitted in ultrarelativistic heavy-ion collisions 
at large relative rapidity has revealed an azimuthal structure that can be interpreted as solely 
due to collective flow [49 ^ I50 [ IST j 152]. This interesting new phenomenon, denoted as triangular 
flow, results from initial state fluctuations and a subsequent hydrodynamic-like evolution. Unlike 
the usual directed flow, this phenomenon has no correlation with the reaction plane and should 
depend weakly on rapidity. Event-by-event hydrodynamics [53j has been a natural framework 
for studying this triangular collective flow but it has been of interest also to investigate these 
correlations in terms of the PHSD model. We have found the third harmonics to increase steadily 
in PHSD with bombarding energy. The coefficient v^ is compatible with zero for ^/sNN > 20 GeV 
in case of the hadronic transport model HSD which does not develop 'ridge-like' correlations. 
In this energy range PHSD gives a positive ^3 due to dominant partonic interactions. 

Different harmonics can be related to each other and in particular, hydrodynamics predicts 
that V4^ oc (^2)^ |44) . In this work it was noted also that V4 is largely generated by an intrinsic 
elliptic flow (at least at high pt) rather than the fourth order moment of the fluid flow. Indeed, 
the ratio ^4/(^2)^ calculated within the PHSD model is approximately constant in the whole 
considered range of ^/sMN but significantly deviates from the ideal fluid estimate of 0.5. In 
contrast, neglecting dynamical quark-gluon degrees-of-freedom in the HSD model, we obtain a 
monotonous growth of this ratio. 

The transverse momentum dependence of the ratio ^4/(1^2)^ at the top RHIC energy has 
given further interesting information (c/. Fig. 4) by comparing the various model results to 
the data from STAR which are interpreted as follows: the STAR data are not compatible with 
ideal hydrodynamics and a finite shear viscosity is mandatory (in viscous hydrodynamics) to 
come closer the experimental ratio observed. The kinetic approaches AMPT and PHSD perform 
better but either overestimate (in AMPT) or slightly underestimate the scattering rate of soft 
particles (in PHSD). 

We recall that our present PHSD calculations employ 'naturally' Glauber type initial state 
fluctuactions which appear compatible with experimental observations up to top RHIC energies. 
On the other hand one might expect that at LHC energies an initial 'glasma' phase might play a 
sizeable role and that the gluon-field fluctuations - of lower scale - could show up in the Fourier 
decomposition of the azimuthal angular distribution. It will be interesting to compare the 
coefficients Vn for high multiplicity pp, p + Pb and Pb + Pb reactions as a function of transverse 



moment um pr (and centrality). 
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